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Abstract. We investigate some properties of the nonequilibrium stationary 
state (NESS) of a one dimensional open system consisting of first and second 
class (type 1 and type 2) particles. The dynamics are totally asymmetric but 
the rates for the different permitted exchanges (10 — > 01, 12 — > 21, and 
2 — > 2) need not be equal. The entrance and exit rates of the different 
species can also be different. We show that for certain classes of rates one can 
compute the currents and phase diagram, or at least obtain some monotonicity 
properties. For other classes one can obtain a matrix representation of the 
NESS; this generalizes previous work in which second class particles can neither 
enter nor leave the system. We analyze a simple example of this type and 
establish the existence of a randomly located shock at which the typical density 
profiles of all three species are discontinuous. 



1. Introduction 

The properties of open systems kept in a nonequilibrium stationary state (NESS) 
through contact with very large (formally infinite) reservoirs at different temper- 
atures and chemical potentials is a central issue in statistical mechanics. Of par- 
ticular interest in these systems are the fluxes of locally conserved quantities, such 
as particle numbers, transported by the system from one reservoir to another. The 
only nontrivial systems for which the NESS is known (more or less) explicitly are 
one dimensional lattice systems with stochastic dynamics corresponding to simple 
exclusion processes or zero range processes [U [2] . The latter case is rather excep- 
tional in that the NESS is described fully by a product measure and thus does 
not exhibit the long range correlations expected to be generic for NESS of realistic 
systems OS]. 

In this note we investigate the particle fluxes and densities in the NESS of the 
open two-species totally asymmetric simple exclusion process (2-TASEP). This is 
a system on a finite one-dimensional lattice of L sites, in which each site is either 
empty (empty sites are also referred to as holes) or is occupied by a type J or a 
type 2 particle. The internal (bulk) dynamics are those in which a type 1 particle 
can exchange with a hole or type 2 particle to its right and a type 2 particle with a 
hole to its right or a type 1 particle to its left; that is, a type I particle at site j can 
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jump to site J + 1 unless that site is occupied by another type 1 particle, and a hole 
at site j can jump to site j — 1 if that site is nonempty. For the moment we allow 
different rates for each type of exchange, setting the time scale by normalizing the 
particle-hole exchange rate to be 1, but in many of the special cases to be considered 
below we will specialize to have equal rates for all bulk exchanges. We will follow 
the convention of [5] and refer to type 1 and type 2 particles as first class and second 
class, respectively, when the bulk exchange rates are all equal, as in the original 
definition of second class particles [6]. Particles enter the system from the left, i.e., 
at site 1; type 1 particles can do so either by filling a hole or by displacing a type 2 
particle, and type 2 particles by filling a hole. Similarly, particles exit at the right 
with the replacement of type 1 particles by either holes or type 2 particles, and of 
type 2 particles by holes. 

The one-dimensional 2-TASEP dynamics have been studied earlier, primarily 
using the so-called matrix ansatz; see [5] for an extensive review of this method 
and its applications. The case in which all exchange rates are equal was considered 
both for the closed system on a ring [7] and for the system on the infinite lattice 
[9] . The case of unequal rates for the system on the ring has also been discussed 
[10[ 111) , with a concentration on the situation in which there is only one type 2 
or "defect" particle. The equal-rates dynamics were also studied, and the NESS 
obtained, for a semipermeable open system, i.e., one in which second class particles 
can neither enter nor leave [HI [T3l [14] ; this restriction of course forces the current 
of second class particles to be zero. In [14] it was shown that the projection of the 
NESS on the second class parties (whose total number was fixed) was given by a 
Gibbs measure in which the interaction was long range (logarithmic) but restricted 
to pairs of neighboring second class particles (that is, pairs separated only by holes 
and first class particles). 

In the present work the restriction to semipermeable systems is dropped. This 
permits us to investigate the dependence of the type 2 particle current on the 
different bulk and boundary rates. Unfortunately, we obtain explicit results only 
under various restrictions on the rates. In Section[5]we discuss results obtained from 
projecting the two-species TASEP onto a one-species TASEP in two different ways, 
a process we call coloring. In Section [3] we use coupling arguments to obtain some 
information on how the currents of the different species depend on the boundary 
rates. In Section [4] we discuss cases in which the system is solvable by a matrix 
ansatz, and in Section l4.1l we describe one such system in which, using the matrix 
ansatz, we establish the occurence of a shock at which all three density profiles are 
discontinuous. 

Table [1] gives our general notation for the rates of the bulk and boundary pro- 
cesses; in some sections below we will adopt simpler notation more suited to the 
restricted rates considered there. We will let Jo, J\, and J2 denote the signed 
currents of the various particle species in the NESS; we always have of course 
Jo + J\ + J2 — 0. Note that always J\ > and Jo < 0, while J2 can have either 
sign. Finally, we remark that, given a specification of the rates, an equivalent sys- 
tem is obtained by interchanging type 1 particles with holes, left with right, u with 
w, and oti with fy, i = 1, . . . , 3. 
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Left (Site 1) 


Bulk 


Right (Site L) 


- 


-> 1 with rate ai 


1 -s- 1 with rate 1 


1 


— > with rate f3\ 


- 


-> 2 with rate ct2 


2 — > 2 with rate w 


1 


— > 2 with rate ^2 


2 - 


-> 1 with rate 


1 2 — > 2 1 with rate v 


2 


— > with rate ^3 



Table 1. Rates for the general open two-species process 



2. Coloring 

In this section we assume that the rates for particle exchanges in the bulk are 
independent of the types of particle involved, that is, we take w = v = 1 in the 
notation of Table [TJ It is well known that for a system with these bulk rates and no 
boundaries, that is, a system on a ring or on the (doubly) infinite lattice, one may 
project the system onto a one-species TASEP system, say consisting of black and 
white particles, in two distinct ways: we always color the first class particles black 
and the holes white, but may color the second class particles either (all) black or 
(all) white; we will refer to the systems arising from these distinct colorings as the 
(1, 2 + 0) system and the (1 + 2, 0) system, respectively. Such coloring arguments 
have been used in the past to understand the two-species TASEP on the ring 
[7] as well as an open semipermeable system (in which the second class particles 
are confined) |14j ; in the latter case, however, they provided only an intuitive 
guide, since the boundary rates for the semipermeable system do not satisfy the 
requirements for consistency with the coloring discussed in the next paragraph. 

To extend fully these colorings to the open system we must assume appropriate 
boundary rates. For the (1, 2 + 0) system we require that first class particles enter 
the system at a rate which does not depend on whether site 1 is occupied by a 
hole or a second class particle, that is, we must take ot\ = 0:3. Similarly for the 
(1 + 2,0) system we must take f3\ = (3%. In this note we consider primarily the 
case in which both colorings are possible (but see Remark l2.2[) . and for notational 
simplicity we will then modify the notation of Table [T] and write the transition rates 
in the following form: 



Left 


Bulk 


Right 


— s- 1 with rate a 
— > 2 with rate 7 
2 — > 1 with rate a 


10 -» 01 with rate 1 
20 -> 02 with rate 1 
12 -» 21 with rate 1 


1 — > with rate j3 

2 — > with rate (3 
1 — > 2 with rate 5 



Table 2. Rates for the case in which both colorings are possible 



When the rates are as in Table [5] we can apply to each of the two colored systems 
the known results for the open single species TASEP [15] , For that system with 
entrance rate a and exit rate (3 the asymptotic (L — ¥ 00) current J of first class 
particles is given by 

(2.1) J(aJ) 



2, if <5,/3 > \ (maximal current region), 

a(l — a), if a < 4 and a < /?, (low density region) 
(3(1 — ft), if < \ and (3 < a (high density region) 
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Similarly, the asymptotic density of particles in the bulk is given by 

2 if <3, /3 > h (maximal current), 

a if a < h and a < f3 (low density), 

linear(a, 1 — a) if a = /3 < | (shock line), 

1 — /3 if /3 < | and /3 < a (high density), 

where "linear(a, 1 — a)" means the profile is linear with densities a on the left and 
1 — a on the right. As is well known, this linear profile occurs due to the formation 
of a shock, which performs an unbiased random walk within the system. Note that 
because the current J is continuous across the shock line a = (3 < 1/2 there is no 
need for a separate entry for this region in (|2.ip . 

Now if the two-species TASEP has rates given by Table [2 then the (1,2 + 0) 
system is a one-species TASEP with entry rate a and exit rate is /? + S, so that the 
current J\ and density p\ in the two-species system are given by J\ = J(a,f3 + S) 
and pi = p(a,/3 + 5), where J and p are obtained from (|2.1[) and (|2.2[) . Similarly, 
the (1 + 2, 0) system has entry rate a + 7 and exit rate (3, so that in the two- 
species TASEP we have Jo = — J(a + 7,/3) and p — 1 — p(a + 7, (3). Finally, 
J% = — J\ — Jo and P2 = 1 — pi — po- Thus the phase diagram of the two-species 
system is completely determined by the phase diagrams of the two colored systems. 
We summarize in Table [3] the currents and in Table [4] the bulk densities in the two- 
species system; the second of these tables has one more row and column than the 
first for the reason mentioned above in regard to (j2.1|) - (|2.2|) . Detailed information 
about the density profile near the boundary of the system is obtainable similarly 
from the results of |15j . 

Remark 2.1. Examination of Tables [3] and [4] shows several interesting features of 
these systems. 

(a) In regions of parameter space in which neither colored subsystem is on its shock 
line the densities are uniform in the bulk. In fact we see from Table 0] that, given 
any p = (/Oo,pi,/02) with J2i Pi = there exist parameters a,/3,^,S which yield 
the pi as the (constant) bulk densities. Recall that for any such p there exists 
a unique steady state p p for the 2-TASEP dynamics on an infinite lattice having 
these as densities. It is then easy to see that in the open system with bulk densities 
p the local state [THl E] at any point a fraction x, < x < 1, of the way through 
the system will be p p ; that is, in the L — > 00 limit the NESS near the site i = [xL\ 
(where [£J denotes the integer part of £) will look like p p . 

(b) When both colored subsystems are in the maximal current region, i.e., when 
Q 1 ,/? > 1/2, second class particles are excluded from the bulk and their current 
vanishes. It follows from the results of [IB] , however, that second class particles are 
present in the system, with a density decreasing as the inverse square root of the 
distance from either boundary. 

(c) When both colored systems are on their respective shock lines, necessarily a = /3 
and 7 = 5 = 0; the system is then a one-species TASEP, again on the shock line. 

(d) When the (1,2 + 0) system is on its shock line and the (1 + 2,0) system is not, 
the latter must be in its high density phase. In that case there is a uniform density 
13 of holes in the system, together with a shock, with random location uniformly 
distributed over the system, separating a region on the left with densities p' Q = ft, 
p[ = a, and p' 2 = 1 — a — (3, from a region on the right with p'$ = /3, p'{ = 1 — a, 




^\ (1,2 + 0) 
(1 + 2,0) \^ 


«,/? + £> \ 
(Maximal current) 


1 r 

(Low density) 


/3 + <5<i,£ + <5<a 
(High density) 


a + 7,/3> ^ 

( A/T 3 "vi in 3 1 r , nT"rprii" I 


Jo — —7, 
4 

J 2 = 


Ji = a(l — a), 

Jo = -4' 
J 2 > 




P < -,0 < a + 7 
(High density) 


*4 

J = -j8(l-j9), 

J 2 < 


Ji = a(l — a), 
J = -0(1-0), 
J 2 sign undetermined 


Ji = OS + S)(1 

Jo = -0(1 - /3), 
J 2 < 


a + j< i,a + 7</3 
(Low density) 


X 


Ji = a(l — a), 
Jo = -(a + 7)(1 - a - 7), 

J 2 > 


Ji = a(l — a) 
Jo = -0(1 - /3) = - Ji 
J 2 -0 



Table [3] Currents for each class of particle. Note that the order of regions in the rows and columns 
is different. Always J 2 = — ( Jo + Ji); the table gives the sign of J 2 , when this is determined. Regions 
marked X cannot be realized, and the region in the lower right corner can be realized only with 
a = f3 and 7 = 6 = 0. 



\ri,2 + o) 

(l + 2,0)\^ 


a,P + S> ~ 
(Maximal current) 


1 a r 

a < -,a < ,3 + £ 
(Low density) 


a = /3 + 5 < - 
(Shock line) 


P + S < \iP + 5 < a 
(High density) 


1 

a + 1,(3 > - 
(Maximal current) 


1 

Pi = 2' 
1 

P2 = 


Pi = a, 

Po= 2> 
1 

P2 = ~ 2 -a 




X 


P < \p < "+7 
(High density) 


1 

Pi = 2' 

Po = /3, 
i 

P 2 =2-^ 


Pi = a, 
Po = P, 
p 2 = 1 — a - (3 


pi = linear(a, 1 — a), 
Po = P, 

P2 = linear(l — a — P, a — (3) 


(H = 1-0-5, 
Po = P, 
P2 = 5 


a + 7 = /3 < i 
(Shock line) 




Pi = a, 
po = linear(l - (3,(3), 
P2 = linear(/3 — a, 1 — a — P) 


pi = linear(a, 1 — a), 
Po = linear(l — a, a), 

P2 = 


X 


a + 7< -,a + -f < (3 
(Low density) 




Pi = a, 
Po = 1 - a - 7, 
P2 = 7 


X 


X 



Table 01 Currents for each class of particle. Note that the order of regions in the rows and columns is different. Regions marked 
X cannot be realized, and the region in which both subsystems are on their respective shock lines can be realized only if a = P 
and 7 = 8 = 0. 



OPEN TWO-SPECIES EXCLUSION PROCESSES 



7 



and p 2 ' = a — (3. The local state a fraction x of the way through the system is then 
the superposition x[i p + (1 — x)fi p . Similar conclusions hold when the (1 + 2,0) 
system is on its shock line and the (1,2 + 0) system is in its low density phase. 

Remark 2.2. One can, of course, also discuss values of the boundary rates which 
respect one but not both of the colorings. Much less can then be said, but there is 
one simple case in which the currents can be completely determined. Suppose that 
ax = "3 = ct, so that the (1, 2 + 0) coloring is possible; then since the exit rate of 
type 1 particles is fii + p2 the current Ji is related to the density pi(L) of type 1 
particles at site L by 

(2.3) J 1 = (0 1 +p 2 )pi(L). 

Now suppose further that P3 =0, so that type 2 particles cannot exit at the right 
boundary and holes can enter there only when a type 1 particle exits. The currents 
across the right boundary (and thus across any bond) must then satisfy 

(2.4) Jo = -ML) = --J±—J u j 2 = -f3 2pi (L) = --^-Ji. 

Pi + Pi Pi + P2 

Note that these results apply for any bulk rate w, although we still require that 
v = 1 in order that the (1, 2 + 0) coloring be possible. 

3. Coupling and monotonicity of currents 

We do not know how to calculate the currents of the various species in our model 
for general rates, but in some cases it is possible to show by means of coupling 
arguments that the currents depend monotonically on the rates. The next two 
theorems give results of this nature. 

Theorem 3.1. Suppose that v = w = 1. Then in the region of parameter space in 
which as > ai and /?3 > p\ : 

(a) The current Ji is monotonically increasing in oti, a 2 , and a^; it is monotoni- 
cally decreasing in p3 and under increase of p\ when fii + f) 2 is constant. 

(b) The (negative) current Jq is monotonically increasing in magnitude as a func- 
tion of Pi) f3 2 , and /?3,' it is monotonically decreasing in magnitude in 03 and under 
increase of ai when ai + a 2 is constant. 

Theorem 3.2. Suppose that v = w = 1. Then: 

(a) If ai = Q!3 then the current J 2 is monotonically increasing in a 2 , in ,$3, and 
under increase of p\ when Pi + fi 2 is constant. 

(b) If Pi — P3 then the current J 2 is monotonically decreasing in f3 2 , in 0:3, and 
under increase of ai when ai + a 2 is constant. 

To develop the proofs of these theorems, given later in this section, we introduce 
the natural order -< 2 -< 1 on the three species of particles, and if r and r' 
are two configurations of our system we write r ^ r' if Tj ^ t[ for i = 1, . . . , L. 
To couple two evolving systems, denoted as primed and unprimed, is to define a 
joint dynamics for the pair (r' '(£),t(£)) in such a way that the induced dynamics 
on each separate system is that of Table Q] with rates v',w', a' it (3[ and v,w,ai,f3 i: 
respectively. We are interested only in order preserving couplings, for which if 
r'(0) t r(0) then r'(t) h r(t) for all t > 0. 

We now define the coupling in the bulk and on the boundaries. 
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Bulk: For this coupling it is necessary that all the bulk rates for the two systems 
be the same: 

(3.1) v = v= w = w=l. 

Under this assumption we define the coupled dynamics by supposing that a Poisson 
alarm clock which rings with rate 1 is associated with each bond, and that when 
the alarm rings any possible exchange across that bond, in either the primed or 
unprimed system, takes place. Equivalently one may describe the dynamics by 
giving rates for the coupled process: if (r/, Tj) ^ ( r l+xi T i+i) then at rate 1, 

( T i n)(T' l+1 ,n +1 ) ->■ (^ +1 ,r i+1 )(r/, n) if t[ > t[ +1 and n y T i+ i, 

(3.2) (t-,t 1 )(t' 1+17 t 1+1 ) -> (t/,t 4+ i)(t- +1 ,t,:) if r[ -< and n h T i+1 , 

( T i t i){t' 1+1 ,t 1+1 ) -> (4i,Ti)(Tj, n+x) if t[ h r' i+1 and r t < r i+1 , 

It is easy to see that no coupling of this general sort will succeed if (|3.1[) does not 
hold. 

Boundaries: As an example we begin by considering the particular case in which 
the state on site L is (t' l ,tl) = (1', 1). We must assign rates to the five allowed 
transitions on site L arising from the exit of the first class particle in one or both of 
the coupled systems, that is, transitions from (l',l) to (0',0), (2',0), (2', 2), (l',0), 
or (l',2); on the other hand, transitions to (0',1), (0',2), or (2',1), which would 
violate the ordering of the two configurations, must not occur. The total rate for 
each transition in each of the two coupled systems must be given by Table [TJ so that 
for example the rate for 1' — > 2', which is the sum of the rates for (1', 1) — > (2', 0) and 
(l',l) ->■ (2', 2) must be/%, the sum of the rates for (1', 1) -> (0',0), (l',l) -> (2',0) 
and (1', 1) — > (l',0), must be /3i, etc. There are thus four equations to be satisfied 
by five rates; it is easy to see that a solution exists with nonnegative rates if and 
only if Pi > /3[ and fi\ 4- /?2 > /3[ + /3' 2 , and that there is then in general one free 
parameter. 

All other boundary transitions may be analyzed similarly; the resulting rates are 
given in Table O In this table x, y, s, and t are parameters which may be chosen 
freely subject to the condition that all the resulting rates be nonnegative, that is, 
to the inequalities 

max{0, 0:2—0:2} — x ^ min{a2, ol'i — oti\, 

< y < min{a2, a' 3 — Qfi}, 
max{0, (3' 2 - /3 2 } < s < min{^, ft - (3[}, 

0<t<mm{/3 2 ,p 3 -p[}. 

The inequalities (|3.3[) will have a solution, and all rates of Table [5] will be nonneg- 
ative, if and only if the parameters of the two systems satisfy 

(3.4) a[ > ax, a[ + a' 2 > «i + a 2 , a' 3 > a 3 , a' 3 > ax, 

(3.5) P'x<Px, P[+P' 2 <Pl+P2, /?3</?3, /3'x<fo. 

In summary: an order preserving coupling of the two systems is possible if (|3.1I) . 
(|3.4[) . and (13.51) are satisfied, and is then given by Table [5j 

Remark 3.3. (a) Note in particular that the marginal rates for boundary tran- 
sitions in either the primed or unprimed system are those of Table [1] separately 
in each boundary configuration. For example, the rate for 1' — > 2' at the right 
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Left 






Right 




XvcLLc 


( u 


0) -¥ 


/ -i / 1 \ 

(I M) 




(■\ 1 


,1) 


-> (0 


UJ 


Pi 


( u 


0) -> 


(I', 2) 




f i ' 


,1) 


->■ (2' 


UJ 


s 


( u 


0) -> 


(I'.O) 


CKi — QIj — X 




1) 


->■ (2' 


2) 


p 2 - s 


( u 


UJ ->■ 


(2 >2) 




I 1 






UJ 


fl fl' 

Pi - Pi - s 


CO' 


0)-> 


(2',0) 


a' — f>o -1- T 


fl' 


1) 


-+(!' 


21 




(2' 


0)-> 


(iM) 


Oil 


(!' 


,2) 


-> (0' 


0) 




(2' 


0)-> 


2) 


y 


(!' 


•2) 


-> (2' 


0) 


t 


(2' 


0)-> 


(i',o) 


«3 - «i - y 


(!' 


,2) 


-> (2' 


2) 




(2' 


0)-> 


(2', 2) 


a 2 -y 


(!' 


2) 


->(!' 


0) 




(2' 


2)-> 


(iM) 


«3 


(!' 


0) 


-> (0' 


0) 


ft 


(2' 


2)-> 


(l',2) 


a 3 — a3 


(!' 


,0) 


-> (2' 


0) 




(!' 


o)-> 


(iM) 


ai 


(2' 


2) 


-> (0' 


0) 




(!' 


0)-> 


(l',2) 




(2' 


,2) 


-> (2' 


0) 




(1' 


2)-> 


(iM) 


a 3 


(2' 


0) 


-> (0' 


0) 


fl 



Table 5. Rates for the coupled process 



boundary is ^ if (t'l> t l) — (l'j 1), as discussed in the example above, and also if 
(r|,r L ) = (l',2)or(r| ) r i ) = (l',0). 

(b) The left-right, particle- hole symmetry described in Section [1] can be extended 
to the coupled system by also interchanging the primed and unprimed systems. 

From now on we will write J- and Ji for the currents in the primed and unprimed 
systems. 

Theorem 3.4. (a) Suppose that the coupling is possible, i.e., that (|3. 1|) (|3.5|) hold, 
and that in addition 

(3.6) p 1 +0 2 =&+&. 
Then the currents of first-class particles satisfy J[ > J%. 
(b) Similarly, Jg < Jo when a± + ai — a[ + 

Proof, (a) By (|3.3|) and (|3.6[) we must have s — (3\ — (3[ = (3' 2 — (3 2 in the coupling, 
so that the rates are zero for the processes (1', 1) — >■ (l',0) and (1', 1) — > (1', 2) at 
the right end of the system. Thus a type 1 particle cannot exit the system unless 
a type V does also, and this implies the result. For (b) one may give a similar 
argument or appeal to the symmetry of Remark l3.3f b). □ 

We can now prove the first monotonicity result given above. 

Proof of Theorem \3.1i We show monotonicity of J\ in a±; all other cases are 
similar. If a.\,a.2,ot3 and /3i,/32,p3 satisfy 013 > ct\ and (3% > f3\, and a[ satisfies 
a 3 > a i > a i, then by defining a' 2 — a 2 , ot' 3 — 0:3, and fi[ = f3i, i = 1,2,3, we 
obtain two systems to which Theorem 13.41 may be applied. □ 

Remark 3.5. We have no evidence that the restrictions in Theorem 13 . 1 1 that 03 > 
ct\ and P3 > P\ in are in fact necessary for the monotonicity; they may well be an 
artifact of the method of proof. 



10 



ARVIND AYYER, JOEL L. LEBOWITZ, AND EUGENE R. SPEER 



We now turn to the proof of Theorem [321 f° r which we introduce a more restric- 
tive coupling with a larger set of excluded pairs. Specifically, the coupling has the 
property that if the initial configuration of the coupled systems contains none of the 
pairs (r-,Ti) = (0',2), (0',1), (2',1), (l',0) or (l',2) then these pairs are also absent 
from configurations (r' '(t), r(t)) for all t > 0. For this coupling the parameters must 
satisfy 

(3.7) u = v! = w = w' = 1, 

(3.8) a-y — u\ = 0:3 = ctg, a' 2 > oil., 

(3.9) #<ft, #+/3 2 =ft+/3 2 , /?3</?3. 

Rates for the coupled process are then given in Table [6] 

Note that in (|3.7|) ()3.9[) we have added conditions to those required for the 
previous coupling (that is, to (|3.1j) . (|3.4|1 . and (|3.5jl ) but also omitted the conditions 
ttj > «i and ft < ft. Note also that fl3.8[) implies that both the primed and 
unprimed systems can be given the (1,2 + 0) coloring discussed in Section [5] and 
hence that J[ — ft ; this also follows from the fact that particles 1 ' and 1 appear in 
the coupled system only in the pair (1', 1). Thus any inequality between J 2 and ft 
implies the opposite inequality between ft and ft. 



Left 


Rate 


Right Rate 


(0',0)->(l',l) 


Oi\ 


(l',l)->(0',0) ft 


(0',0) -)■ (2', 2) 


a 2 


(l',l)->(2',0) £ 2 -ft 


(0',0)->(2',0) 


a 2 - a 2 


(l',l)->(2',2) ft 


(2',0)-*(l',l) 




(2',2)^(0',0) /3' 3 


(2',0)->(2',2) 




(2',2)->(2',0) ft -ft 


(2', 2) -*■(!',!) 


a 3 


(2',0)-+(0',0) /3' 3 



Table 6. Rates for the coupled process giving monotonicity of ft 



Theorem 3.6. Suppose that (f3T7j) - ([33)) fto/d. T/ien: 
fft // ft = ft for i = 1,2,3 tfien ft > ft . 
(b) If a.2 = a' 2 then J 2 < ft- 

Proof. Under the hypothesis of (a) no (2', 0) pair can enter the system at the 
right boundary, while under (b) no (2', 0) pair can enter the system at the left 
boundary. □ 

Proof of Theorem VS. 2c (a) is obtained form Theorem 13.61 as in the proof of Theo- 
rem 13.11 and then (b) follows from the symmetry of Remark 13.31^ . □ 

One can check easily that these results are confirmed by those of Table [3J when 
both colorings are possible; in that case they imply that ft should be increasing in 
/? and 7 and decreasing in a and S. 

4. Matrix solutions 

In this section we observe that for various special choices of the boundary rates 
of Table Q] the probabilities of configurations of the open two-species TASEP can 
be computed using a variant of the usual matrix product ansatz. Specifically, this 
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means that for these rates there exist matrices Xo, X%, and X2 and vectors \V) 
and {W\ such that the probability of the configuration t\ . . . tl, where tj is 0, 1, or 
2, is given by 

(4.1) Mi( T l, •■•,tl) = Z^WLin, . . . ,t l ); 
here the weight wl is given by 

(4.2) w l {t x ,...,t l ) = (W\X T1 X T2 ---X TL \V), 
and the normalizing factor Zl by 

(4.3) Z L = w L {Ti,n,...,T L ) = (W\(X + X 1 +X 2 ) L \V). 

Ti,...,T L 

Such matrix solutions exist for two different families of boundary rates, which we 
call permeable and semipermeable. We can impose either type of boundary condition 
at either end of the system, leading to four distinct families of open systems. 

We now describe the rates (summarized in Table [7]) and the necessary algebraic 
properties of the matrices Xi and vectors \V), (W\ for these systems. 

Bulk: The exchange rates in the bulk are those of Table [IJ although we shall on 
occasion specialize to take v and/or w equal to 1. For the weights ()4.1[) to be 
stationary the matrices Xq, Xi, X2 should satisfy 

_ 1 1 

(4.4) X1X0 = Xi + Xo, X1X2 — —X2, X 2 Xq = —X2, 

v w 
the usual algebraic rules for the 2-TASEP [TD1HT]. 

Semipermeable boundaries: A semipermeable boundary is one through which 
type 2 particles cannot pass, that is, a semipermeable left boundary has a 2 = ct3 = 
and a semipermeable right boundary P2 = /?3 = 0. A semipermeable boundary is 
thus characterized by a single rate, which we will write as a = ot\ for a left boundary 
and /3 = /?i for a right boundary. The vectors (W\ for the left boundary and \V) 
for the right boundary should satisfy 

(4.5) Semipermeable left boundary: (VK|A = — (Wl, 

a 

(4.6) Semipermeable right boundary: -^l|V) = ttI^O- 

P 

The open system with two semipermeable boundaries and with w = v = 1 was 
studied in [IB]. [T2]. and [H]. 

Permeable boundaries: At a permeable boundary type 2 particles can enter and 
leave the system, but the rates are constrained by 

(4.7) Permeable left boundary: ol\ + a 2 — w, 

(4.8) Permeable right boundary: /?i + /?2 = v. 

We will choose a convenient parametrization of such rates which will be discussed 
further below. At the left boundary we write a 3 — a for the rate of exit of type 2 
particles; there is a second free parameter, z, which plays the role of a fugacity for 
type 2 particles, and a.2 = waz, a\ = w(l — az). At the right boundary we write 
similarly ^3 = 6, $2 = vbz, and /3i — v(l — /3z). If both boundaries are permeable 
then the parameter z must be the same on the left and the right, and to insure 
nonnegative rates these parameters must satisfy 

(4.9) az < 1, bz < 1. 
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The vectors vectors \V) and (W\ used to construct the steady state should satisfy 

(4.10) Permeable left boundary: (W\X = -{Wl (W\X 2 = z(W\, 

w 

(4.11) Permeable right boundary: X t \V) = ~\V), X 2 \V) = z\V). 





Left 
SP 


p 


Bulk 


Right 

SP P 


-> 1 


Q 


w(l — az) 


10 


-> 01 


1 


1 ->• 


/3 w(l-6z) 







waz 


20 


-^02 


w 


1 -» 2 


vbz 


2 -4 1 





a 


12 


-►21 


V 


2^0 


b 


(W| 


(w a \ 


(W w \ 




IV) 


IK) IK) 



Table 7. Rates for permeable (P) and semipermeable (SP) 
boundary conditions. The last line shows the boundary vectors 
(W\ and \V) to be used in (|37Tj) . 



It is a straightforward exercise, using standard arguments, to verify that if there 
exist matrices and vectors Xq, X\, X 2 , \V), and (W\ satisfying (|4.4I) . either (|4.5[) or 
(|4.10| . and either (|4.6[) or (|4.11j) . then (|4.ip defines an invariant measure for the cor- 
responding dynamics. We establish this existence below. One striking consequence 
is that if the system has one or more permeable boundaries then the parameters a 
and/or b, which control the entry and exit of type 2 particles, do not appear in the 
algebraic conditions (|4.4[) and (I4.10P and/or (|4. 1 1 [) . so that the steady state of the 
system is independent of these parameters. 

To establish the existence, recall [15] that there exist matrices Xq and X\, and 
vectors (W a \ and \Vp) satisfying 

(4.12) X 1 X a = X 1 +X , 

(4.13) (W a \X = -(W a \, X l \V f} ) = hvp) 

a p 

(4.14) (W a \V )= "f a + (3>l. 

a + p — 1 

We now define 

(4.15) X 2 = z W + V ~ 1 \V v )(W w \; 

wv 

then the matrices Xq, X\, X 2 satisfy (|4.4I) (in fact, this requires only X 2 = 
X\V V )(W W \ for some non-zero A PHIE])- Note that for a system with two semiper- 
meable boundaries we are now introducing an additional parameter z (which is 
again a fugacity, as discussed below). For semipermeable boundary conditions on 
the left (respectively right) we take (W\ = (W a \ (respectively \V) = \Vp), thus 
satisfying (I4.5[) (respectively (|4.6|1 ). For permeable boundary conditions on the left 
(right) we take (W\ = (W w \ (\V) = \V V )); (|3TU)) (respectively (|4TTT|l ) will then be 
satisfied due to the choice of overall constant in the definition (|4.15l) of X 2 . This 
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completes the construction in the case when all the matrix elements (|4.2|) finite 
which, in view of (|4.14j) . is certainly true in the region 

(4.16) a + f3>l, a + v>l, w + (3 > 1, and w + v > 1. 

To extend the discussion to all parameter values we must distinguish two cases. 

Case 1: One or both boundaries permeable. In this case all weights (|4.2[) have 
the form QP{\/a, 1/(3, 1/w, l/iO, with P a polynomial having positive coefficients 
and Q a common factor depending only on the type of system under consideration: 
Q = (w + v — in the case P/P (that is, two permeable boundaries), Q = 
(w + P - I)- 1 in the case SP/P, and Q = (a + v - l)" 1 in the case P/SP. The 
factor Q cancels in the probabilities (|4. 1 [) and the resulting probabilities are well 
defined and positive for all positive values of the parameters; the algebraic identities 
among them which express stationarity hold by analytic continuation from the 
region (|4.16l) . 

Case 2: Both boundaries semipermeable. In this case the number of type 2 
particles in the system is conserved by the dynamics, and the NESS decomposes 
into ergodic components corresponding to the L+l possible values of N2 . Within 
each component an argument as in Case 1 allows one to extend the NESS to all 
parameter values; however, the factor Q referred to there depends on N2, with 
Q = Q = (a + f3 - l)- 1 if N 2 = and Q = Q* = (w + v - l)(a + v- l)~ l (w + 
(3 — l) -1 if A^ > 0. One may also consider grand canonical ensembles in which 
N2 fluctuates. The algebraic conditions (|4.4p - (|4.6j) do not determine the relative 
weights of configurations with different values of N2, but the particular choice 
of representation of the algebra discussed above does, and thus defines a grand 
canonical ensemble in the region (I4.16|) via (14. (I4.3[) . In general the inequality 
of Qo and Q* is an obstruction to extending this ensemble to all values of the 
parameters — negative probabilities can arise. This difficulty can be avoided, and 
a grand canonical NESS defined, by excluding the component with N% = from 
the superposition. Alternatively, a grand canonical NESS including all components 
can be defined when Qo = Q*, that is, when w = a or v = j3. 

Remark 4.1. (a) The cases of semipermeable and permeable boundary conditions 
are not mutually exclusive; on the left, the case a = of permeable conditions 
corresponds to the case a — w of semipermeable, with a similar common case on 
the right. 

(b) The system with permeable boundaries on both ends was studied, using other 
methods, in |17j . The authors took (in our language) a = b = vw/(l + vwz) and 
did not observe that the NESS is in fact independent of a and b. 

As noted above, the case of two semipermeable boundaries with w = v = 1 was 
studied, in the canonical ensemble with fixed N2, in [13] and [14], and it is easy to 
see that several features discussed there are still present for arbitrary v and w. The 
current J 2 is always zero. The operator X2 has rank one (see (|4.15[) ) which implies 
that a type 2 particle at site i decouples sites to the left of i from those to the right, 
except that the total number of type 2 particles in these two subsystems must be 
N2 — 1. As a consequence [14] the marginal measure on configurations of type 2 
particles is an equilibrium measure for a certain Hamiltonian with nearest-particle 
interactions. In any of the grand canonical ensembles discussed under Case 2 above 
a second class particle completely decouples the system; here we use the fact that 
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the grand canonical ensembles for the left and right subsystems are well defined 
(see the last remark under Case 2). 

Let us consider then the case in which there are permeable boundary conditions 
at one or both ends of the system. As noted above the NESS is then independent 
of the parameters a and/or b, and so is determined by v, w, z, and a (respectively 
ft) if semipermeable boundary conditions are used at the left (respectively right) 
boundary. In particular, we may achieve this steady state with a = or b = 0, so 
that certainly we still have J2 = 0. That is, all models with the rates of Table [7] 
have zero current of type 2 particles. 

In all of the grand canonical ensembles discussed above the parameter z plays 
the role of a fugacity for the type 2 particles. To see this, suppose that r is a 
configuration with N2 type 2 particles and note from (|4.5p and (|4.6j) that the weight 
assigned to r by the algebra above will contain a factor z N2 . Now let wl,n 2 ( t ) 
denote the weight assigned to r by the algebra (|4.4p - (|4.6[) of two semipermeable 
boundaries with a = w and j3 — v (for v — w — 1 this is the weight assigned to r 
in the algebra of [14]). Then we see that 

L 

Hl{t) = Zl 1 z N2 w l , N2 (t). 

N 2 =0 

That is, fiL is just the "grand canonical" ensemble, with fugacity z for the type 2 
particles, for the semipermeable system with entry and exit rates w and v, respec- 
tively, in the sense that this term is often used when discussing matrix ansatz mod- 
els. It then follows, for example, that the type 2 particles are in a grand-canonical 
version of the equilibrium ensemble mentioned above. 

4.1. A simple special case. We close this section with the consideration of a 
special case in which various quantities can be computed rather easily: we take 
permeable boundary conditions on the left and semipermeable on the right, and 
specialize to v = 1 (there is a corresponding analysis with left and right reversed 
and w = 1). The algebra then reduces to 



X\Xq — X\ + Xq, X1X2 — X2, X2X0 — — X2, 



(4.17) 



Xi\V)=^\V), {W\X = -{W\, (W\X 2 =z(W\, 
p w 



and is realized by taking \V) — \Vp), (W\ — (W w \, and X 2 = z\V\){W w \ If we now 
define new matrices X[ and X' 2 by X[ = X\ and X' = Xq + X2 then these satisfy 

[ 1 -4- in 7 

(4.18) X[X' Q =X' 1+ X' Q , X' 1 \V) = -\V), (W\X' Q = ——{W\. 

This is the algebra of an open one-species TASEP system [TS] with entry rate 
a = w/(l + wz) and exit rate j3, which implies that if we color the type 1 particles 
black and the holes and type 2 particles white, as for the (1, 2 + 0) system discussed 
in Section [51 then the marginal of the NESS on the black/ white system is precisely 
the NESS of this 1-TASEP system. (The situation here is somewhat different from 
that of Section [21 however; the correspondence here is on the level of measures 
rather than dynamics, since the rate at which a type 1 particle enters the system 
depends on whether there is a hole or a type 2 particle on site 1.) 
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Many properties of the 2-TASEP may now be obtained. The current and density 
profile of type 1 particles are given by (|2.1|) and (|2.2|) : 



(4.19) J 1 = J(a,/3), Pl =p(a,j3). 

The current of type 2 particles is zero since the right boundary conditions are 
semipermeable, and the density profile of type 2 particles, 

(4.20) p 2 (x) = lim Prob(r l = 2), 

may be obtained from the fact that X 2 is a rank one matrix: 

(W\(X +Xx +X 2 y- 1 X 2 (Xo+X 1 +X 2 ) L - i \V) 



Prob(r 4 = 2) = 
(4.21) = z- i_1 L ~ l 



(W\(X Q + X 1 + X 2 ) L \V) 



and the known asymptotics [15] of Z^' 5 = (W 7 |(Xo + X\) L \Vs). 

As an example we discuss the case a — (3 < 1/2; for convenience we will eliminate 
z and /3 and use a and w as basic variables to describe the system in this case. We 
will need the asymptotic formulas [15] 

L(l - 2 7 ) 2 



if 7 = <5< 1/2, 
if 7 < S and 7 < 1/2. 



(<5-7)a L (l-a) L+1 
In this case all density profiles are linear: 

(4.23) ^ = linear (r-i, i = 0,l,2, 
where 

, , (1 — 2a) (w — a) a(l — 2a + aw) 

(4.24) n = a, r 2 = ^ n ; , r = 

w(l — a) w(l — a) 

(4.25) ri = 1 - a, = 0, r' = a. 

For pi this follows from (I4.19p . for p 2 from (|4.20|) and (I4.2ip . using z — (w — a)/(wa) 
and (|4.22l) . and then for po from YliPi = 1- 

We now show that these linear profiles arise, as usual, from the occurrence of a 
shock in the typical (not averaged) density profiles: a randomly located point at 
which (in general) all three typical profiles are discontinuous. In this case the shock 
location is marked by the position of the rightmost type 2 particle in the system; 
the probability that this particle is located at site i is, from (|4.22|) . 

M9fi v (WKgo + x 1 + x 2 y-^x 2 {x + xtf-^v) _ z^z^ i 
1 ' (wkxo+x^iv) z z2' a ~v 

so that this particle (and hence, as we will see below, the shock) is uniformly located 
throughout the system. When conditioned on the location i of this particle the 
NESS decomposes into the product of the NESS's of two independent subsystems, 
and by finding the densities in these subsystems we exhibit the discontinuity in the 
overall (conditioned) density profile that is the characteristic of a shock: 
Left subsystem: On sites 1 . . . , i—1 the NESS is that of a semipermeable/permeable 
2-TASEP like the full system, but with exit rate /3 = 1. One may again introduce a 
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(1,2 + 0) coloring for this system for which the type 1 particles form an open system 
with entrance rate a and exit rate 1; the colored system is in its low density phase 
[15] and thus has (up to boundary effects) a uniform density of type 1 particles: 
p-y = a — n . The density profile p2 can then be calculated from (|4.20[) and (|4.21|) 
with the replacement (3 — ¥ 1; the result is a constant profile with value rg. The 
profile po is then also constant with po = ro- 

Right subsystem: On sites i+1 . . . ,L the NESS is that of an open one-component 
TASEP with entry rate w and exit rate /3 = a. Since a < w this system is in its 
high density phase; density profiles are constant with p\ = 1 — a and po = a. 

Note that if w = 1 then r$ = r' Q , so that the density of holes is constant through 
out the system and the densities of first and second class particles have compensat- 
ing discontinuities at the shock position. Situations of this kind were encountered 
in Section [5] (see Remark l2.1f d)) and in the study of the "fat shock" in the system 
with v = w = 1 and two semipermeable boundaries |14j . For general w all three 
densities are discontinuous at the shock positions, which so far as we know is a new 
phenomenon. 

The shock is stationary in the sense that it has no net drift velocity, although its 
position does of course fluctuate. Let us write Ji(p; v, w) for the current of species i 
in a uniform steady state of an infinite system with densities p — (p±, P2, P3), where 
the exchange rates in the bulk are 1, v, and w, as in Table [TJ The condition for a 
stationary shock as described above is then that 

(4.27) Ji(r;l,w) = Ji(r';l,w) fori = 0,1, 2. 

It would be useful to have analytic expressions for the currents Ji(p;v,w) so that 
these relations could be checked; such expressions were obtained for the case v = 
w = 1 in [7], essentially by a coloring argument. However, for any given p, v, and 
w one may compute Ji(p;v,w) with good accuracy by Monte Carlo simulation of 
the uniform system on a large ring; we have done so and thus verified that (|4.27[) 
holds for a variety of choices of a and w. 

Remark 4.2. A simple mean field calculation of the currents in a uniform two 
species system yields Ji(p;v,w) = pi(pa + VP2), J2(p',v,w) = p2(wpo ~ v Pi), and 
Jq{P]v,w) = po(—pi — WP2). These mean field values are in fact correct when 
v = w = 1, as is well known, but the Monte Carlo simulations referred to above 
show that they are not correct in general. 
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